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ABSTRACT 

9 

The cross section for photoneutron production in Be at energies of 
1 . 692, 1 . 720, and 1 . 78 Mev has been investigated using radioisotope 
gamma-ray sources. Neutron counting was by manganous sulfate bath 
and a long counter. Gamma-ray counting was by a scintillation spec- 
trometer consisting of a Nal(Tl) scintillation crystal and a 256-channel 
pulse ^-height analyzer. The following results were obtained: 

a j = 1 . 35 ± 0. 08 mb 
= o. 763 ± 0. 022 

a 1. 692 



— - ■ = 0. 617 ± 0. 013 

a 1. 692 

A theoretical treatment of the P-*-S transition near threshold using 
square-well potentials is also presented. 

I. INTRODUCTION 

ft 

The photoneutron cross section in beryllium is of great interest for 
several reasons: (1) beryllium has the lowest threshold for this reaction 

of any element; (2) beryllium is used in many photoneutron sources where 
a supply of monoenergetic neutrons is desired; (3) beryllium is often used 
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as a reflector in nuclear reactors; (4) the photodisintegration cross sec- 
tion may be used as a test for models of the beryllium nucleus. The pho- 
toneutron reaction was discovered in 1934 by Chadwick and Goldhaber 
(Ch34), who observed the reaction in deuterium, and in the same year 
Szilard and Chalmers (Sz34) observed the reaction in beryllium. Early 
measurements of the cross section for photoneutron production in beryl- 
lium were attempted by several experimenters, with widely varying re- 
sults. The cross section was calculated theoretically by Guth and Mullin 
(Gu49) in 1949., and their results have been fairly well accepted until very 
recently. The most recent theory is essentially similar to that of Guth 
and Mullin, but substitutes a potential of the Saxon type for the square well 
originally used, resulting in a considerably better fit with experimentally 
determined points. 

This paper will present a brief discussion of photoneutron reactions 
in general, some discussion of the theory, a discription of several new 
experiments measuring the cross section at different energies, and con- 
clusions drawn from these experiments. 

II. GENERAL 

Neutrons were first discovered in 1930, having been produced by the 
(a, n) reaction. The (y, n) reaction was then discovered in 1934, when 
neutrons were observed to result from the irradiation of beryllium and 
deuterium by high-energy gamma rays. Further investigation showed 
that these photoneutrons were emitted with a well-defined energy spectrum, 
in contrast with those resulting from alpha-particle bombardment, which 
have a widely distributed energy spectrum. In many applications, this 
discrete energy characteristic is very useful. Unfortunately, however, 
there is an attendant disadvantage. The cross section for the production 
of photoneutrons is very small, so that in order to produce a reasonable 
neutron flux a gamma-ray source of very high intensity must be employed. 
Unlike the alpha particles, which are stopped within the source, the gamma- 
ray intensity is attenuated only negligibly, thus producing high gamma-ray 
backgrounds and making handling conditions difficult. 
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(y, n) reactions in general involve the absorption of a gamma photon 
into the nucleus, forming an excited nucleus which then decays with the 
emission of a neutron. Thus 

= E^ — (binding energy). 

From this it becomes apparent that if the gamma-ray energy is well- 
defined, the neutron energy will also be well-defined. The best source 
of monoenergetic gamma rays is the radiation from radioisotopes. Un- 
fortunately, the useful gamma rays (that is, those having a long half-life 
and high yield) are limited to energies of less than 3 Mev. Although the 
(y, n) reaction occurs in many elements, only deuterium and beryllium 
have thresholds below 3 Mev for this reaction. By the use of various 
radioisotopes these two elements have been employed in sources for the 
production of monoenergetic neutrons with energies from 25 to 1000 kev. 
Photoneutrons can also be produced by the use of br ems s trahlung from 
targets bombarded by high-energy electrons from an accelerator, but 
these neutrons are not monoenergetic due to the continuous energy spec- 
trum of the bremss trahlung. 

III. THE (y, n) REACTION IN BERYLLIUM 

Photoneutrons are produced from beryllium by the following reaction 
for gamma-ray energies less than about 2. 5 Mev: 

9 8 

hv + Be -►Be + n — 1 . 666 Mev 

g 

In this reaction the Be is left in its ground state and decays with a half- 
—1 5 

life of ~10 second as follows: 

Be^ -► 2 He^ + 0. 09 Mev 

It would appear that the reaction 

9 4 

hv + Be -► 2 He + n — 1 . 576 Mev 

is favored energetically, since it has a lower threshold. In 1952 Alburger 
(A152) looked for this reaction, using 1.631-Mev gamma rays from Ne , 
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and observed no neutrons which could be attributed to this reaction. On 

the basis of this experiment he concluded that the cross section for the 

—3 0 2 

reaction could not be greater than 10 cm . 

At gamma-ray energies higher than about 2. 5 Mev many other re- 
actions are possible. These have been summarized by Paolini in a re- 
cent paper (Pa60) and will not be repeated here. 

9 8 

The total cross section for the Be (y, n)Be reaction has been meas- 
ured experimentally at several energies using radioisotope sources. 

Some of the more important of these are listed in Table I. 



TABLE I 

Measurements of Photoneutron Cross Sections for 
Beryllium Near Threshold 



Gamma 

source 


Gamma energy 
(Mev) 


Cross section 
(mb) 


Reference 


Sb 124 


1. 692 


1 . 0 mb 


Note 1 






1.9 


Note 2 






1. 262 ± 0. 069 


(Gi59) 






1.35 ± 0. 08 


This paper 


Bi 206 


1. 720 


1. 03 ± 0. 07 


This paper 


Al 28 


1. 78 


0 . 83 ± 0.06 


This paper 


Y 88 


1 . 85 


0. 654 ± 0. 031 


(Gi 59) 






0. 55 


Note 3 


o 144 
Pr 


2. 185 


0. 39 


(Ha 53) 


MsTh 


2. 61 


0. 50 


(Al 50) + (Bi50) 






0. 377 ± 0. 02 


(Ed56) 


Na 24 


2. 754 


0. 70 


(Ru48) 






0. 674 


(Sn50) 



Notes: 

1. Unpublished. Found by A. H. Snell about 1946. 

2. Obtained from (Ru48) by application of branching ratio. 

3. Unpublished. Found by E. Segre in 1949. 
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In 1956 Connors and Miller (Co56) explored the cross section near 
threshold using bremsstrahlung from an accelerator and reported a peak 
at 1. 70 ± 0. 05 Mev, with the cross section > 200 times the average cross 
section in this energy region, and with a peak width < 50 kev. A more re- 

cent experiment by Walter, Shea, and Miller (Wa60) fails to show such a 
peak. In this experiment they also used bremsstrahlung from an electron 
accelerator, using a thin gold target, and varied the maximum electron 
energy in 5-kev steps. Their results agree with the cross section as 
measured with radioisotope sources. They analyzed the bremsstrahlung 
spectrum using the Sauter-Fano theory, finding that this gave a better fit 
than did the earlier Bethe -Heitler theory. They attribute the peak reported 
in the earlier experiment to an incorrect analysis of the bremsstrahlung 
yield from the thick gold target used at that time. 

A recent experiment at M.I. T. by Paolini (Pa60) determined the over- 
all photoneutron spectrum from about 3 to 18 Mev. This is shown in Fig. 1. 
In this figure the peak just above threshold as determined with monoener- 
getic sources has been added to Paolini f s results. 

IV. THEORY 

In 1949 Guth and Mullin (Gu49) calculated a theoretical cross section 

9 

for photoneutron production in Be , based on a two-body model consisting 

g 

of a single neutron loosely bound to a Be core, and using square-well 
potentials. In their model the ground state is taken as a P state * Elec- 
tric dipole transitions then lead to an ^^/z state near threshold and to D 
and states at somewhat higher energies. The cross section for a 

magnetic dipole transition, which would lead to a ^ ^ /z s ^ ate near thresh- 
old, was calculated to be very small. Their expression for the cross sec- 
tion for the P — S transition contains four adjustable parameters namely, 
the radius and depth of the initial and final wells. Making the assumption 
that the initial and final wells have the same radius of action and choosing 
the initial well depth to give the correct binding energy, reduces the num- 
ber of adjustable parameters to two. For a given radius the depth of the 
final well is then chosen to result in a fit to the experimental cross section 
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near threshold. In their calculations they used a radius of 5. 0 f (1 fermi 
— 1 3 

= 10 cm) and a final well depth of 3 . 0 Mev. The threshold (binding 
energy) at that time was thought to be 1. 63 Mev. The total cross section 
calculated by Guth and Mullin is shown in Fig. 2. 

An early theoretical treatment by Bergmann (Be5l) was similar to 
that of Guth and Mullin, but the numerical results were different. Some 
of the data included in his paper indicated the effect of varying the pa- 
rameters . 

It appeared to us that the use of more recent experimental data, to- 
gether with a judicious choice of parameters, might produce a better fit 
to the experimental cross section. On this basis, new square-well cal- 
culations following the treatment of Guth and Mullin were made. The de- 
tailed development of this theory is included in Appendix A of this report. 
Calculations were made using several values for the radius of the initial 
and final wells, and the initial well depth in each case was chosen to give 
the correct value of 1. 666 Mev for the binding energy of the last neutron. 
The depth of the final well in each case was chosen for two different con- 
ditions: (1) correct energy dependence of the cross section, based on the 

experimental ratio of the cross section at 1. 78 Mev to that at 1. 692 Mev; 
and (2) correct value of the cross section at 1. 692 Mev as experimentally 
determined. The results of these calculations are shown in Figs. 3 
through 5, which also illustrate the effect of changing the various pa- 
rameters. It is felt that, because of the magnitude of the parameters 
used to fit the experimental cross section, the square-well assumption is 
not entirely valid. The qualitative fit, however, is quite good, especially 
considering the simplicity of the model. It is interesting to note that the 
sharp peak just above threshold demands that the final state be a lightly 
bound state ( ~50 kev binding). 

Goldman, Francis, and Guth (Go60) have recently calculated the 
theoretical cross section near threshold for the P — S transition using a 
single -particle shell model of the Saxon type: 



cr x 1 0 28 (cm 2 ) 
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9 

Fig. 3. Be (y,n) cross section curves illustrating the effect of 
varying r 0 . Vq j a> is chosen to fit the ratio of the cross sections 
at 1. 692 and 1. 78 Mev. Note that the absolute magnitudes have been 
multiplied by reduction factors. 



(mb) 
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Fig. 4. For three values of r Q , Vq has been chosen to fit the 
absolute value of the Be^^.n) cross sfecxion at 1. 692 Mev. 



Q ~ ps (mb) 
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9 

Fig. 5. The Be (y,n ) cross section curve with r and V.. . > chosen 
to fit the absolute magnitude at 1. 692 Mev and the ratio of tne'cross 
sections at 1. 692 and 1. 78 Mev. 
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V (r ) 



-V 0 

~ r-“R 

1 + exp — — — 



with V 0 determined from the binding energy of the last neutron, and with 
the final well based on scattering data. Using R = 3. Of for both the in- 
itial and final states, with a = 0. 6 f for the initial state and a = 1. 2f for 
the final state, they have obtained a cross section curve which appears to 
fit the experimental data in both energy dependence and absolute magnitude. 



V. EXPERIMENTAL PROCEDURE 

124 

A. Sb Gamma-Ray Source 

9 

It was decided that the total photoneutron cross section of Be for the 

124 

1. 692-Mev gamma rays from Sb would be measured. This was the 

closest gamma-ray energy to the threshold which was readily available in 

sufficient intensity and with a long enough half-life to permit convenient 

measurement. Also, this energy appears to be quite near the top of the 

124 

peak in the cross section near threshold. Furthermore, Sb sources 

had been used in previous experiments (see Table I) which had resulted 

in cross section measurements differing by almost a factor of 2, and it 

was desired to clear up this long-standing discrepancy. 

A source consisting of a hollow sphere of 1 / 1 6 in. aluminum packed 

with granulated antimony to a density of about 2. 7 g/cm , and having a 

1.3 6-cm outside diameter, was irradiated in the Livermore Pool Type 

Reactor (LPTR) to a total antimony activity of 22 5 millicuries and allowed 

to decay for a period of several months to eliminate any short-lived activity. 

This gamma-ray source was then held within a hollow beryllium sphere 

made up of nesting hemispheres, the whole assembly of beryllium having 

an outside diameter of 8.2 cm and a thickness of 0. 645cm. A sketch of 

the assembly is shown in Fig. 6. The particular source geometry used 

was chosen to: (1) give spherical symmetry; (2) minimize the correction 

1 24 

necessary for the finite si^e of the Sb source; and (3) reduce the ab- 
sorption of gamma rays and the moderation of neutrons in the beryllium 
shell to a reasonable amount. 
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Fig. 6. Beryllium photoneutron source. MUL-8456 
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B. Neutron Counting 

In order to measure the neutron flux, a long counter was constructed 
following the pattern of the Harwell 4 model (A155). The effective face of 
the counter was determined by using Po-Be, Pu-Be, Po-Li, Sb-Be, and 
mock fission sources at various distances, and plotting (count rate) 
vs distance from the actual face of the counter. The use of the long 
counter for absolute neutron flux measurements proved unsatisfactory, 
however, because of the differe nt geometries of the various sources and 
an uncertainty as to the variation in efficiency of the counter over such a 
broad range of neutron energies. 

As a better means of counting the total neutron flux, a manganous 

sulfate bath was prepared (Ge59). This method of determining neutron 

flux has the advantage that essentially all the neutrons are moderated and 

captured by the manganese, so that the resulting activation is independent 

of the neutron energy spectrum. A tank 44 in. in diameter and 29 in. deep 

was filled with a nearly saturated solution of manganous sulfate, using 

175 gallons of distilled water in which was dissolved 575 pounds of MnSO^. 

A Victoreen 1B85 Geiger-Muller tube was attached to a probe of aluminum 

137 

tubing and waterproofed for dip counting. A Cs button source was fixed 

in a standard geometry for periodic checks of the G-M tube counting rate. 

An NBS Ra-Be source was used to calibrate the MnSO. bath. It was sus- 

4 

pended in the center of the tank in a waterproof holder made of aluminum 

5 6 

and the bath was activated to saturation. Since the half-life of Mn is 

2. 58 hours, any activation over a period of a day or longer essentially 

saturates the bath. The G-M tube was then placed in the center of the 

bath and a series of counts taken over a period of several half-lives of 
56 

the Mn to ensure that the proper activity was being counted. From the 
results of this series of counts, the count rate at the time of removal of 
the source from the bath was computed. From this initial counting rate 
and the total flux of the NBS source, the total flux of any' other neutron 
source can be determined by immersing it in the bath and activating the 
bath to saturation. The ratio of the fluxes is then the same as the ratio 
of the initial counting rates. The initial counting rate obtained with the 
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NBS source has proved to be reproducible to within 1% in subsequent 
measurements. This procedure was followed in determining the total 
neutron flux produced by the antimony-beryllium source assembly. 

C. Strength of Gamma-Ray Source 

In order to compute the photoneutron cross section, it was necessary 
that the total flux of gamma-ray photons with an energy of 1. 692 Mev be 
determined. This was accomplished by the use of a thallium-activated 
sodium-iodide scintillation crystal and a Z 56-channel pulse-height an- 
alyzer. Computations of the gamma-ray intensity were based on un- 
published curves for intrinsic crystal efficiency and photofraction fur- 
nished by Harry West of this laboratory. It is estimated that the ac- 
curacy of this method is better than 5%. The scintillation counting tech- 
nique is felt to be the best method of determining the strength of the 
gamma-ray source in this experiment, since it is independent of the de- 
cay scheme of the source, and provides an automatic correction for self- 
absorption and Compton scattering in the source. In using other methods 
the strength of a particular gamma line must be calculated from a know- 
ledge of the decay scheme involved, and in many instances the decay 
scheme is not known with great accuracy. With the scintillator and 
pulse-height analyzer combination, only those gamma rays that have a 
specific energy are counted, and further, any gamma rays absorbed in 
the source are not counted, just as they are not seen by the beryllium 
surrounding the gamma-ray source. In addition, any gamma rays which 
are Compton scattered are not counted under the photopeak. This is de* 
sirable, in the case of antimony, since almost all (99. 5%) of those photons 
which have undergone Compton scattering are below the threshold for 
photodisintegration. Thus the gamma-ray source looks the same to the 
scintillator as it does to the beryllium, and corrections for absorption 
and Compton scattering are applied automatically, giving a more ac- 
curate evaluation of the gamma-ray flux available to produce photo- 
neutrons . 
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D. Computation of Cross Section 

Once the gamma-ray flux (at the required energy) and the neutron 
flux are known, the total cross section for photoneutron production at 
the particular energy involved can be computed from 



where 



a = 



N A 
n Be 



N y No p 



Be 



Ar 



N 



n 



Be 

N 

V 

N 0 



p Be 

Ar 



cross section in cm 

neutron flux in neutrons/sec 

9 

atomic weight of Be 
gamma-ray flux in photons/sec 
Avogadro* s number 

9 3 

density of Be 7 in g/cm 
9 

thickness of Be shell in cm. 



E. Comparison Technique 

It was desired to explore the photoneutron cross section near thresh- 
old in as much detail as possible using radioisotopes, but the gamma-ray 
sources available in this energy region are limited to a few which have 
short half-lives and/or low intensities. With a half-life less than several 
hours or a total neutron flux less than about 10 n/sec, the use of the MnSO 
bath becomes subject to large uncertainties. Therefore, a new experi- 
mental arrangement providing for the following considerations was de- 
sirable : 

1. The measurement should be independent of half-life. This could 
be accomplished by a simultaneous count of gamma rays and neutrons. 

2. It should be possible to count a weak neutron flux. This could be 
done with a long counter. 

3. The effect of room scattered neutrons should be minimized. This 
could be done by placing the source close to the long counter. 

4 

4. Since there are about 10 gamma photons for each neutron prod- 
uced, and since the scintillation counter is more efficient than the long 
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counter, precautions must be taken to prevent excessive gamma count- 
ing rates, with attendant dead-time losses. This requirement could be 
satisfied by locating the scintillator at a distance from the source, and 
by using a collimator at the crystal as necessary. 

The arrangement finally devised is shown in Fig. 7. The simultane- 
ous count of neutrons and gammas results in an n/y ratio which is pro- 
portional to the cross section. When several gamma sources of differ- 
ent energies are used successively in the same geometry, a relative cross 
section curve is readily obtained. This can then be converted to an abso- 
lute cross section curve if the correct aboslute value is known at one of 
the energies. 

In computing the n/y ratios, corrections are made for the variation 
in intrinsic efficiency and photofraction of the Nal scintillator. No cor- 
rections were made for variations in efficiency of the long counter. This 
counter was of the Harwell 4 type which Allen found to be flat down to 
Sb-Be photoneutron energy (A155). Even for the typical 10-15% droop of 
the efficiency of a long counter, our efficiency uncertainty would not be 
greater than about 2%. The long counter was biased to eliminate y count- 
ing. 

This comparison technique was used to make cross-section measure- 

2 8 

ments at 1. 78 Mev with A1 as the gamma-ray source, and at 1. 72 Mev 

206 

with Bi . A sphere of solid aluminum the same size as the antimony 

source was activated in the LPTR and placed in the center of the beryl- 

2 8 

lium shell. It was necessary to work rapidly, since the half-life of A1 
is only 2.3 minutes. This was accomplished by removing the source from 
the reactor by a pneumatic tube delivery system, so that the total elapsed 
time after removal from the reactor flux until the count of the Al-Be as- 
sembly was started was less than one half-life. A lead collimator having 
a 3/4 -in. aperture was used at the scintillation crystal to cut down the 
counting rate for the gamma rays while still maintaining a good neutron 
counting rate. Several successive runs were made as the source decayed, 
and those runs having good counting rates without excessive dead-time 
corrections were used in computing the n/y ratio. The n/y ratio for 
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antimony in the same geometry was measured and the cross section at 
1 . 78 Mev computed by 

* n / V A . 78 

~(n/y) 1.692 ai - 692 “ CTl - 78 ' 

It should be noted that the ratio of the cross sections is considerably 
more accurate than the absolute value. 

For the measurement at 1. 72 Mev a foil of high-purity lead was bom- 
barded in the 60-inch cyclotron at Berkeley to produce the reaction 
pb206 ^ 2n)Bi^^. The procedure was then the same as that used in 
the case of the aluminum source, except that the bismuth source was 
considerably weaker, so the collimator was not used at the scintillation 
crystal. The ratio n/y for the bismuth source was somewhat more dif- 
ficult to determine than for either the aluminum or the antimony, since 
some Bi was present from the Pb (d, 3n) Bi reaction. A rather 
complex gamma-ray spectrum resulted from this contamination, and it 
was necessary to separate out its effect on the basis of pulse height and 
half-life. 

There are several other radioisotopes with gamma-ray energies 
above threshold and short half-lives which could be used in this compar- 
ison technique. We know of no others, however, which have gamma-ray 

2 8 

energies between threshold and the 1. 78 Mev from the Al . Use of 
higher energies than this would probably demand an individual calibration 
of the long counter efficiency curve. 



VI. EXPERIMENTAL RESULTS 



In measuring the photoneutron cross section at 1. 692 Mev, using 



gamma rays from Sb 



124 



the following results were obtained: 



Strength of NBS Ra-Be source E-1277 as of 9 December 1954: 



Q 0 = (6. 09 ± 3%) x 10 b n/sec. 

Strength as of 1 7 September 1959 considering growth of Po 



210 



Q = Q„ 



1 . 12 — 0 . 12 exp 



( — t/ 28)] 



= (6.21 ± .3%) X 10 n/sec. 
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(dN/dt) 0 for NBS source in MnSO^ bath = 1292 c/min 

(dN/dt) 0 for Sb-Be source assembly = 11.05 c/min 

N from Sb-Be = }*’ ° 5 X 6.21 X 10 6 = 5.31 X 1 0 4 n/sec 
n i l p a. 

for 1. 692-Mev line of Sb as of 3 September 1959 = 5. 61 X 10 
as of 17 September 1959 = 5.61X10^ exp ( — 0.6 93 X 14/60. 9) 

= 4. 78 X 10^ photons/ sec 
The following corrections and uncertainties apply: 

Correction to cross section Correction (%) Uncertainty(%) 



Gamma-ray source strength 5. 0 

Neutron counting statistics ~ 2.0 

Finite source size —0.3 <0.1 

Attenuation of gamma ray in Be shell +2.5 <0. 1 

Effect of 2. 09-Mev gamma ray —4. 4 1.0 

Neutron escape from MnSO^ bath —0.9 <0.1 

Absolute strength of NBS source E — 1277 3. 0 

Total correction —3. 1 

Overall uncertainty 6. 2 



The corrections for finite source size, attenuation of the gamma ray 
in the Be shell, effect of the 2. 09-Mev gamma-ray, and neutron escape 
from the MnSO^ bath are described in Appendix B. 

The cross section is now computed as follows: 



a 



a 



N 

n 



N N 
Y 




Ar 



(1 - 0. 031) 



5. 32 X IQ 4 X 9. QZ X Q. 969 

4. 78 X 10 8 X 6.02 5 X 10 23 X 1. 85 X 0. 645 
1. 35 X 10~ 27 cm 2 
1 . 3 5 ± 0.0 8 mb . 



In determining the cross section at 1 . 78 Mev the comparison tech- 

2 8 

nique described previously was used, with A1 as the source of gamma 



rays, and the following results were obtained: 

n/y for Sb-Be = 0. 1179 
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..206 



n/y for Al-Be = 0.0728 

0. 0728 v , , c _ A Q , __ u 
a 1.78 0.1179 X 1 ' 35 0. 83 mb . 

The ratio of the cross sections at 1. 78 Mev and 1. 692 Mev is consider- 
ed to be much more accurate than the absolute value of the cross section, 
having an estimated uncertainty of ± 2%. This gives an overall uncertainty 
at 1. 78 Mev of 6. 6%, so that 

0*1 - 0.83± 0.06 mb . 

In the case of the cross section at 1 . 72 Mev, the strength of the Bi‘ 
gamma-ray source was rather low, so the collimator was removed from 
the sodium iodode crystal to obtain better gamma counting statistics. 

This change in geometry changed the n/y ratio for the Sb-Be assembly. 

The following results were obtained: 

n/y for Sb-Be = 0. 0421 
n/y for Bi-Be = 0. 0309 . 

From these figures a preliminary cross section at 1.72 Mev was 
computed as 



a 



1. 72 



~ 0.0309 
0. 0421 



X 1. 35 = 0. 991 mb 



It was necessary to make an additional correction to this prelimi- 

.205 

nary value, since there was an appreciable amount of Bi present in 

20 5 

the source. Bi has strong gamma rays at energies of 1. 766 and 1. 777 

2 06 

Mev. These, together with the 1. 72-Mev gamma ray from the Bi 3 
resulted in a broad photopeak made up of one peak at 1. 72 Mev and an- 
other, approximately one-fourth as intense, at about 1. 77 Mev. Even 
if the ratio of the two peaks were 1:3 or 1:5 rather than the estimated 
ratio of 1:4, the resulting cross section would not change by more than 
about 1 . 3%. The assumption of a cross section of 0. 85 mb as a reason- 
able value of 1. 77 Mev leads to 



so that 



0. 8 a 1 72 + 0. 2 (0. 85) = 0. 991 mb; 



a j = 1.03 mb . 



Because of the complicated spectrum, an uncertainty of ± 3% was 
estimated for the comparison procedure in this case. This leads to 
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an overall uncertainty of about 7% for the cross section at 1.72 Mev, 
with the result that 

a i 7Z0 = (!♦ 03 ± 0. 07) mb . 

In regard to the absolute value of the cross section, it should be 
mentioned that the antimony source was also counted by the Hazards 
Control Group using a Lands werk Condenser R-meter. The total 
strength of the source was measured as 47. 5 me as of 1 1 August 1959. 

The strength as of 17 September 1959 was then 31.2 me. If we assume 
that the 1. 692-Mev gamma ray occurs in 50. 6% of the disintegrations 

g 

(Gi59), the strength of this gamma ray was 5. 84 X 10 photons/sec. This 
must be corrected for scattering within the source, since scattered 
gammas are counted by the R-meter but are below the photoneutron 
threshold. The correction is about 6. 5%, so that 

N = 5.46X10^ photons/ sec . 

Y 

Subsituting this value in the cross section formula we obtain 

0* = 1 . 1 8 ± 0.07 mb . 

As stated previously, it is felt that this method of measurement is not 
as reliable as that using the scintillation crystal. However, the difference 
in the values obtained by the two methods indicates the presence of some 
systematic error, possibly in the branching ratio. 

In an attempt to improve the accuracy of the original cross section 
measurement at 1. 692 Mev, particularly the neutron counting statistics, 
a stronger antimony source has been prepared and the neutron flux 
measured. This source was also counted by the Hazards Control Group. 
Their measurement of the stronger source leads to the same cross sec- 
tion, within 1%, as did their measurement of the earlier source. The 
estimated uncertainty for the neutron counting statistics as shown in the 
table of corrections to the cross sections is based on the measurement 
of this stronger source. Even though the cross section based on scin- 
tillation spectrometry is considered more accurate than that based on 
R-meter measurements, the latter have proved consistent in successive 
measurements of an antimony source, within the known accuracy of the 
half-life . 
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VII. CONCLUSIONS 

The value obtained in the present experiment for the photoneutron 
cross section at 1. 692 Mev is in reasonable agreement with the recent 
accurate determination by Gibbons et al. (Gi59). This agreement re- 
moves the long-standing large uncertainty in the value of the cross 
section at 1 . 692 Mev. 

The energy dependence of the cross section near threshold is now 
well defined by the combination of our results at 1. 692, 1. 720, and 1. 78 
Mev with the data of Gibbons et al. (Gi59) and Walter et al. (Wa60). 

The use of a square-well potential in the theory of Guth and Mullin 
does not give a fit to the data unless the radius is made unrealistically 
small (1. 87 f). The new calculation by Goldman et al. (Go60) using a 
Saxon-Woods potential with reasonable parameters results in a good fit 
to the experimental data. 
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APPENDIX A 

THEORETICAL DEVELOPMENT USING 
SQUARE-WELL POTENTIALS 



This development follows that of Guth and Mullin (Gu48) but is 
presented in more detail. We considered it necessary to verify the 
formulae of Guth and Mullin because of the different numerical re- 
sults obtained by Bergmann (Be51). 

The general formula for the electric dipole photoneutron cross 
section is 



a 



i 



|2 

M (E) i 

z | 



( 1 ) 



in which 



M (E) = Ze 
z ' m 



I 




ik dr. 

T i 



(2) 



[i/m is the ratio of the reduced mass of the system to the mass of the 
residual nucleus, z = r cos Q is the electric dipole operator, and 
E = hv — € , where € is the binding energy of the neutron to the residual 



nucleus . 

Using a ground P state, the selection rule Al = ± 1 for electric 
dipole transitions permits P -> S and P-*D transitions. Only the P S 
transition will be treated here, since it is the one which is important 
near threshold. 

9 

The expression for the electric dipole disintegration of Be then 
becomes 



a = 




■£ z dr 



(3) 



where the summation extends over all non- vanishing matrix elements. 
Assuming the spin interaction to be small and of a central (non-tensor) 
type, the eigenstates f, i are described by the quantum numbers E 
(energy), j, rm, and 1. The eigenfunctions are 
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^ (1, j = 1 + 2 . m j) = 



R i..i (r> 

(21 + 1) 1/2 



< 1 + k* m j )1/2 



Y. 1 ( 9, <b ) a 

1. m. -r r n 
J- z 



+ (l + I-m)‘/ 2 Y + l( ft< ,)p n 

J 2 



(4a) 



1 R 1 j< r > 

^ (!. j = 1 -T. m ;) = — 175 

J (21 + 1 ) ' 



(l + i-tn ) 1/2 Y j _1 <9.*>a 

J J 2 



; fi * 



-(l + I + m ) l/2 Y 1 <9,*)P n 

J J 2 



(4b) 



The matrix element J z i|^ df of Eq. (3) is different from zero only if 

Al = ± 1, Aj = 0, ± 1, and Am. = 0. a and (3 are spin operators with 

$ $ J n ^ $ 

the property that a a = 6 8 = 1 and a 8 = 3 a =0. 

n n 'll m n 'll n n 

For a P 3/2 i^bial state, ny. can be ± 3/2 or ± 1 /z with equal prob- 
ability. Since the final S-^/z state can on ly have rm = ± l/2, a factor of 

l/2 in the final cross section results from the selection rule Am. = 0. 

3 

Furthermore, Eq. (4b) would apply to a P state, so that only Eq. (4a) 
need be considered in the case of a ^ initial state. 

Considering m. = + l/z, the initial eigenfunction is 

^1 3 /z^ T ^ ( 1 /? 

= ' ’ ( 3 ;i/2 1 (2) ' Y l,0^>a n + Y 1|]L (0^)P n j. (5a) 



4; (1, 3/2, + 1/2) = 



For m. = — 1/2, the initial eigenfunction is 

R, . /,(r) 



+ (1, 3/2, -1/2) --—fir- ’Y 1| _ 1 (e,«.)a n+ (2) 1/2 Y l o ( 0 ^)P n 



(5b) 



For rm = + l/2, the final eigenfunction is 

f 

43 (0, 1/2, +1/2) = R 0j l/2 (r) ' Y 0 0 (e,f) a n 

For m. = — l/z, the final eigenfunction is 
J 

■HO, ‘/2, -1/2) = R 0 ,l/2 (r) | Y 0,0 (e - > P„ 

r * { 

The expression \ i|^ z 4 ^ dT can now be written as 

1/2 rrr , o * \ v \ _2 



(2/3) 



(5c) 



(5d) 



* 

Jz x o,o 



Y„ „ \ r cos Q (R, Y, n ] r'" sin 9 dr d0 d^ 



“1,3/2 1,0 
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Separating variables gives 

(2/3) l/2 J R* l/2 R t> j/2 r 3 dr JJ Y* „ Y^ „ cos 8 sin 0 d 0 d* . 

The normalized spherical harmonics are 

Y - fiJl '/ 2 Y = MV / 2 
0,0 ',471 J 1,0 \4tt/ 



which results in 



( 2 ) 

4tt 



1/2 



J R 0, 1/2 R l, 3/2 r3 dr j'j' COs2 esiMd8if 



Integrating the angle part gives 

,!/2rv* 



( -|l J ^0, 1/2 "T, 3/2 



r 3 dr 



R 0 , 1/2 R 1 

Considering the factor of l/2 mentioned previously, the expression 



now become s 



| J iJj* z i)j. d T 


2 






2 

R 


i 2 = i 


R 


L 3 _ 


| ps 


9 


ps 



where 



r * 3 

R = \ R n i /-> R i i /-> r dr 

ps J 0,1/2 1,3/2 



This makes the expression for the cross section 

2 
R 



3 2 

1 28tr & 

° = ~~7Z9 — V 



ps 



By using hv = and h = 2 tt]^, this becomes 

i »’2 



a = 



64tt 2 e 2 
729 fic 



/w R 



ps 



The only remaining problem is the evaluation of R 
spherical square-well potential, for the initial state 

V = -V, 



ps 



(7) 

Using the 



1, 3/2 



r<_ r 0 



0 



r> r 0 



and the radial wave functions are 
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R 1 , 3/2 



= B 1 + ar e ^ (r-r o > 
°1 2 e 
r 



r > r . 

— o 



where 



and 



with 



a = 



p = 



% « 

% (V 



1/2 



1, 3/2 



+ E) 



1/2 



= reduced mass of neutron - 

1 dR 



8M 



Matching the logarithmic derivative ^ at r = gives 

A 1 

= — sin |3r . 



pr Q cot pr 0 = 1 + (1 + ar 0 )(£) . 

r * 00 22 

Normalizing by using \ (R^ 2>Jz^ r dr = 1 gives 

(2 + ar Q ) (| ) 

+ (1 + ar Q ) (|) -1 



^ f 1 + 

2(3 2 



2 

sin (3r 



o U 



= 1 



(P*o> 



For the final S state, assuming the same radius of action for the 
initial and final states 



V = - V 



0 , 1/2 



r < r 
— 0 



V = 0 

The radial wave functions for this state are then 
R 



r > r 

— o 



■ 0 , 1/2 



‘0, 1/2 



A o 


j 0 (Y*) 


r 


A 


sin yr 




0 


V* 




B o 


$in Tk(r - r 0 ) + 6 0 J 


r 



— 0 



> r 

— o 



V = 



■jt <Vi/a + E > 



11/2 



where 
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and 






1/2 



It will be noted that the wave function for r > r is that for a free 

— o 

particle with phase shift 6 • 

o 

Matching the logarithmic derivative at r = r leads to 



A o 


sin yr Q 


= /- B sin 
k o o 


A o 


cos vr 
Y 0 


= B cos 6 • 

0 0 



The normalization in this case is not so .straightforward as for the 
initial state. It is accomplished as follows: 



I ' 



(R 



0 , 1/2 
r 



,2 2 , , 
) r dr = 1 



I 



r) 



r 2 dr + B 2 

0 



i 



sin j^k(r — ^ ) + 6 J 



kr 



( 8 ) 



r dr = 1. 



The first integral is readily integrated: 



0 J o 



sin ~y r 
Y r 



2 , A 

r 2 dr = -4, 

2Y 



|yr o - cos Y r 0 sin Y r 0 



3 



To integrate the second integral, assume for the moment a small energy 
spread such that dk = (k^ — k^) << 1 and write the integral 

2 r”f Sin |>l (r ~ r . )t6 o] sin [ k 2 (r - r . ) + 6 j l r 2 dr . 

• J l 



k l r 



k 2 r 



For the upper limit, substitute g and rewrite as 

.2 „g 



lim r s 

\ sin Tk, (r — r ) + 6 *j sin fk (r — r ) + 6 | dr . 

g-^qo J r L 1 0 oJ 1-2 o oj 



k l k 2 g ' 

Now substitute x = r - r^ , giving 



B 



lim 



g_>0 ° 



n S 

\ sin k^x+6 0 ) sin (k^ x + ) dx 
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This can now be integrated to give 




Now, since (k^ — k^) ^ terms resulting from the left side of Eq. (8) 

are negligible except 



B 



lim 



g ^°° 



sin (k 1 - k 2 ) g' 

2 < k i-v _ 



This can be rewritten in terms of the Dirac delta function as follows 
(see L. I. Schiff, Quantum Mechanics, McGraw-Hill, 1959, p. 51): 




For normalization per unit dk interval then 

,2 



B 






or 



B q = k(2/ir) 



1/2 



For use in Eq. (7), however, j R 



psl 



must be normalized per unit 



energy interval. This is accomplished by multiplying B q by 



so that 




B 



/dk\l/2 _Ajtk 2^/2 

°W V 



In order to determine R it is now only necessary to subsitute the 

ps 

normalized wave functions into Eq. (6) and integrate. This is a lengthy 
but straightforward integration, resulting in 
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? r~ 7 ~1 ^ / 2 

B 1 a [ (2/tt) (jik/j< )J 


2 + 2 (p 2 /a 2 ) (1 + ar Q ) - r 2 (y Z - (3 2 ) 


n 2 2 ~2 .1/2 
r Q (k + y cot yr Q ) 


(V 2 - p 2 ) 2 



ar (2 + ar ) + k 2 r 2 
o o_ o_ 

/ 2 j. i2\ 2 

(a + k ) 



3 + ar o + (k 2 /a 2 ) (1 + ar Q ) 

/ 2 , , 2.2 

(a + k ) 



2 + [<Y 2 -p 2 ) A 2 ] <1 + ar„ ) 



77^77 



vr cot vr 
Y o Y o 



j 



APPENDIX B 



CORRECTIONS TO THE CROSS SECTION 

It is necessary to make several corrections to the cross section 

computed from the basic formula 

N A 
n Be 

CT " F1I 



y 0 P Be Ar 



These will be described in the following paragraphs. 



1. Finite source size. We follow the derivation given by Gibbons 
et al. (Gi59) „ 




-33- 



UCRL-5996 



Assuming a homogeneous spherical source of strength Q and radius 

3Q 

a, the strength per unit volume is »• 

4na 

If Nt = number of atoms/cm in the Be shell, then 

,2 



, » T . R 3Q 

dn = Nta — j 

P 4iTa 



dV 



gives the number of neutrons per second due to source element dV. The 
total neutron flux produced is r^f where n^ = Nt<jQ is the result if the 
source were concentrated in a point at the center of the shell, and f is 
the finite source size correction given by: 



R 2 fdV _ 

V J 2 



3R 

. 3 

47ra 



III 



r sin 0 dr d 0 d 0 
2 2 

R + r — 2Rr cos 0 



This integrates to 
3R 
2a 3 



f = 



Ra - 



R 2 —a 2 



In 



R + a^ 
R — a> 



Since a<< R, the logarithm may be expanded: 

x— 1 , 1 [x — 1? , 1 |X — 1 



In x = 2 



x + 1 



Letting 



x = 




(x> 0) 



R - a 



x + 1 



In x = 



2a 

IT 



7~ 


l / < 


i\ 2 1 , 


^a\ 4 1 


1 


+ 3 \B 


0 + 5^ 


vRy + ' ’ * 



Subsituting in the_equation for f: 

4 



f = 



3 



2a 

f = 1 + 



2 2,2 

3 a + TT 



+ - 



R 



5 



5R 

For extreme accuracy this should be multiplied by the average value of 
sec a, but this correction is negligible when a << R. 

Z. Attenuation of gamma ray in Be shell. 
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3. Effect of the 2.09-Mev gamma ray. 

In the case of the*Sb-Be source assembly, the neutron flux measured 
in the manganous sulfate bath includes not only the neutrons produced by 
the 1. 692-Mev gamma ray, but also those produced by the 2.09-Mev 
gamma ray. A correction for this effect must be made in computing the 
cross section at 1 . 692 Mev. The ratio of the 2. 09-Mev gamma ray to 
the 1. 69-Mev gamma ray was measured as 0. 148 ± 0. 008 by Nal count- 
ing. A cross section of 0.4 mb was assumed as a reasonable value at 
2. 09 Mev. The correction to be applied to the cross section is then 

0. 148 X 0. 4 = 0. 06 mb. 

This was divided by the computed cross section to give a percentage 
correction. Our correction is the same as that applied by Gibbons et al. 
(Gi59). 

4. Neutron escape from MnSO^ bath. 

Neutron escape from the MnSO^ bath was significant only in the case 
of the NBS standard source. This correction was computed by dividing 
the outer boundary of the bath into zones, computing the neutron flux in 
each surface zone, and weighting by the solid angle subtended at the center 
of the bath. The fraction f of neutrons escaping in each zone was com- 
puted by the formula 

f = -f = 1. 5e _0 * lr , (Ge59) 

0 o 

where r is the distance from the source to the edge of the tank in 
centimeters . 
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This report was prepared as an account of Government sponsored work. 
Neither the United States, nor the Commission, nor any person acting on 
behalf of the Commission: 

A. Makes any warranty or representation, expressed or implied, with 
respect to the accuracy, completeness, or usefulness of the information con- 
tained in this report, or that the use of any information, apparatus, method, 
or process disclosed in this report may not infringe privately owned rights; or 

B. Assumes any liabilities with respect to the use of, or for damages 
resulting from the use of any information, apparatus, method or process dis- 
closed in this report. 

As used in the above, M person acting on behalf of the Commission n 
includes any employee or contractor of the commission, or employee of such 
contractor, to the extent that such employee or contractor of the Commission, 
or employee of such contractor prepares, disseminates, or provides access 
to, any information pursuant to his employment or contract with the Commis- 
sion, or his employment with such contractor. 



